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STEERING LAW FOR PARALLEI MOUNTED 
DOUBLE -GIMBALED CONTROL MOMENT GYROS 

i. INTRODUCTION 


There are several reasons why some space vehicles require angular 
momentum exchange devices in addition to (or maybe instead of) reaction con- 
trol systems: Fine pointing, long mission duration, and lack of contamination. 

FI le pointing requires a restoring torque which is continuously variable over 
several orders of magnitude; whereas, reaction Jets are on-off devices. Lang 
mission duration can be handled by momentum exchange devices since they use 
electrical energy which can be replenished by solar cells, but reaction lets use 
mass expulsion and the fuel consumption over long periods becomes prohibitive. 
Momentum exchange devices will not contaminate the immediate space vehicle 
environment, in contrast to reaction jets. 

Three types of momentum exchange devices are commonly used: reaction 
wheels ( RW ' s ) , single-gimbaled control moment gyros (SG CMG's), and 
doublc-glmbalod control moment gyros (IXI CMG's). RW's achieve momentum 
exchange by accelerating (or decelerating) about a fixed axis. Therefore, they 
have one degree of freedom ami their energy content varies greatly, requiring 
relatively large average power. Wheel imbalance can generate disturbances 
at all frequencies. The momentum steering laws are simple, therefore only 
satellites or smaller space veh'cles use RW' s. 

SG CMG's have constant momentum magnitudes which can lx: positioned 
anywhere in planes perpendicular to the individual gimbal axes. For a system 
of SG CMG's, singularities (where no control torque can be generated in some 
direction) are a proolem, dictating either a complicated steering law or 
the inefficient use of tin: angular momentum. 

1X1 CMG's have iwo degrees of freedom each which allows maximum use 
of the constant momentum magnitudes and the singularities are easily avoided. 
Maximum momentum usage coupled with minimum software requirements also 
mandates unlimited angular freedom of the gimbal axes (no gimbal stops and 
the associated problems). When these IXi CMG's are mounted with their outer 
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gimbai „xes parallel and when the steering law 1 can accommodate any number 
of CMG's, a generally applicable momentum exchange system can be devised 
which always has a spherical momentum envelope. The fact that the number of 
IXJ CMG's is also freely selectable allows a few standard sizu CMG's to satisfy 
any angular momentum requirement; however, IXJ CMG's will be used mainly 
on large space vehicles (Skylab, Space Shuttle). 

Figure 1 shows the inner and outur gimhul angles o ^ and /i . of the 
i-th CMG. With this definition, the total angular momentum of n CMG's is: 
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where i = 1,2, ... , n and II. is the momentum magnitude of the i-th £MG. 
The CMG momentum change is 


IF, * I A|d ♦ | Bl/J , ( 2) 



Figure 1. 1X1 CMG gimbai angles. 


1. The term steering law means here the sel of equations net essary to gen- 
ei ate a set of gimbai angle rate commands, which then result n the desired 
control torque on the space vehicle. This is achieved by steering the momentum 
vector of each CMG according to the commands generated by the steering law. 
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Thu torque exc*rtu«l on the CMG .system is 
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and it is convenient to define the control torque command T as 
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is the desired momenlt. m change. 


4 


The mounting of all outer glmbal :ixom parallel to the Home axis (dis- 
cussed in detail in Section II) makes the control torque (generated along this axis 
only a function ol the inner ((I mlju 1 angles and their rutes (eqs. 1-13). This 
fact allows decomposition of the total three-dimensional control problem Into a 
linear one for the inner rate commands and a planar one for the outer rate 
commands (Section IV). Since there are more than three degrees of freedom 
for the total system, an infinity of choices exists for the gimbal rat commands 
while satisfying the commanded control torque. The pseudo-Inver' method is 
selected to give the basic sot of gimbal rate commands. A two stei approach is 
user'. The inner rate commands are calculated lirst, and they are hen con- 
sidered known quantities In the calculation of the outer gimbal rati commands, 
both times using the pseudo-inverse method. 

However, since the pseudo-inverse does not prevent sing* lar states, u 
set of non -torque -producing inner and outer glmbal rate comma' ds are added to 

the basic set (Section III), such that all singularities (besides aturation) are 
avoided. 


Proportional glmbal rate command limiting is appll' I to all commands 
in the case any one of the gimbal rate commands exceeds no respective torquor 
capabilities. This ensures that only the magnitude of the generated torque is 
reduced, but the direction of the torque vector is not altered. 

This report Is the result of simplifications mentioned In Reference 1, 
but it was possible to carry the simplification so far as to eliminate the need 
for the variable weiguts of tin; pseudo-inverse which wore used to invert the 
nonsquare lAl and fit) matrices. Therefore this report presents actually a 
different approach to the steering law problem. 2 


II. PARALLEL MOUNTING ARRANGEMENT 


The pro|>o.scd mounting, with all outer glmbal axes parallel, is shown 
in Figure 2. Although the outer glmbal axes are shown eolinear, this is not a 
requirement. This mounting arrangement lias many advantages. The mounting 
interfaces can l>c identical, i.e. , mounting brackets and hardware, cable 
harnesses, etc. There is no need to individually identify the IXJ CMC’ s, and 
the onboard computer can assign an arbitrary label lo any CMG, which could Iks 
different from one computation cycle to the next. This simplifies the steering 
law and the redundancy management. The parallel mounting of 1 he outer gimbal 
axes in conjunction with a steering law that accepts any number of IXi CMG'h 
also makes failure accommodation a built-in feature. On the other hand, if 


2. The steering law was called "control law" in Reference I. 



Figure 2. Parallel mounting arrangement. 



increasing momentum requirements during the design of a vehicle (and the 
moments of inertia always lend to Increase) demand it, an additional CMd ^un 
‘ n added with minimum impact on hardware and software. The parullel mount- 
„ of 1X1 CMli' s also makes visuali/allon of system 0 |»er»itio.. *xccedingiv 
simple (especially when compared will) the momentum envelope of skewed 
mounted single gimbal CMC’ s) . 

The only disadvantage of the parallel mounting lies in the fact that a 
large momentum demand along the direction of the outer gimbal axes reduces 
the maximum response in a plane perpendicular to the outer axes, since a small 
movement of the momentum vector demands a large outer gimbal rate. This 
possible problem can lie eliminated by mounting the outer gimlxil axes in the 
direction of the minimum momentum requirement (usually the minimum prin- 
cipal moment of inertia axis). 

III. DESIRABLE GIMBAL ANGLE DISTRIBUTIONS 


The n IXJ CMOs have 2n degrees of freedom. Three are needed Jo 
satisfy the torque command; the excess of 2n-:i degrees of freedom are utilized 
to achieve a desirable gimlml angle distriixition. Before one can decide on a 
desirable distriixition, the characteristics of a IKI CMG have to Iks considered. 

The inner gimbal rate needed to produce a given torque |>crpondiculnr to 
the inner gimbal axis is independent of the ii net* and the outer gimlxil angles. 
However, the outer gimbal rate needed to produce the same torque pcr|H?ndlcular 
to the outer gimlxil axis is Inversely proportional to the cosine of the inner 
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gimkil angle. Therefore It is desirable tu keep the cosines of the inner gimlxil 
angles high; l.e. , it is desirable to minimize the muximum inner gir.ibnl angle, 
which, in turn, reduces the outer rate requirements. The maximum inner 
gimbal angle is minimized when all inner gimbal angles arc equal (also com- 
pare equation ( l) ) to the inner gimbal reference angle 

ir - .In " 1 I n«, VI "i I • ('“) 

n \ r v ‘ ‘ y r 1 j 


The simplicity of equation ( IK) is the direct result of the parallel mounting. 

The situation is not as clear-cut with respect to the desirable distribution 
for the outer gimbal angles. However, for IXJ CMGs, a singular condition inside 
the total momenta > envelope can only occur when some of the vectors (at least 
one) are antiparallel and the others parallel to the resultant. Maintaining 
adequate and more or less equal spacing between the vectors will eliminate the 
possibility of a singularity. In this report, a distribution is used which maxi- 
mizes the minimum separation between any two outer gimbals. The angles Ad 
between adjacent vectors are then equal, as shown in Figure :i for three, four, 
and five DG CMGs (for two DG CMGs, the only excess degree of freeoom is 
used for the equalizatb" of the inner gimbal angles). Each of the vectors is 
halfway between its neighbors (the extreme vectors being obvious exceptions). 
The resultant momentum in the V 2 - V 3 plane then bisects the angle between 
the two extreme vectors for ‘.he case where all momentum magnitudes are 
equal and the inner gimbal angles are equal. 



Figure :i. Outer gimbal angle distribution. 


The calculation of the value ior involves the solution of u transcen- 
dental equation, which in addition has a different form depending on whether the 
number of CMG' s is odd or even. These difficulties can ije avoided by the 
generation of a set of reference angles which lie '.a if way between the 

M 1 

adjacent vectors: 
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and 
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where k. is generated by an ordering scheme such that the value of k, is the 

index of the largest d. . ko is f'.e value of the next smaller one, etc. Details 

of the logic needed can be seen from the flowcharts in Appendix A. In the 
flowcharts, the possibility of one or more failed CMG' s (H ■ 0) is also taken 

into account by replacing n by n-m , where now n is the original number of 
CMG' s and m Is the number of failed CMG’ s. 


IV. INNER GIMBAL RATE COMMANDS 


The change of the V, angular momentum component is not a function 
of the outer gimbal rates ( sec equations 1 through 11) , 
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After a summation over j is done, an inner gimbal rate command, o , of the 
form 
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( 23 ) 
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will result In T„ , » T,,, , If It is assumed that the actual and the commanded 

Cj 1 Cl 

glmbal rates are equal. 

After a summation over i is done, the torque equivalence can l>e shown 
by 
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The effect of the inner gimbal rales on the \ 2 and V t momentum components 
is treated later. 

Simultaneously with satisfying the torque command, it is desirable to 

reduce the maximum inner glmbal angle as discussed in the previous section. 

The inner glmbal angles can be moved toward equality without a V, torque 

l>eing generated, as shown in Appendix It, by an inner gimbal rate addition, 

<i , of the form 
Ai 
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Ai 


K a (•* - “'j) 


( 25 ) 


where 


O' = 



II eo ( II . eo .) 

yr 1 1 
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and K Is a constant gain ( l/K Is the time constant) . The angle <» closely 
A A 

approximates the Inner glmbal reference angle n , as shown in Appendix B. 


The combined Inner glmbal rate command is then* 


(i 




T Cl' K A (< ' -°|» 


( 27 ) 


Rate limiting will lx.- discussed after the outer gimbul rate commands have been 
established. 


V. OUTER GIMBAt RATE COMMANDS 


Again the total rate command is split into a rate addition for vector 
distribution and a rate to satisfy the torque command. 

T!;'* preliminary outer glmbal rate additions are 



k a <"m 



(2K) 


with the exception of the two outermost vectors which will have zero preliminary 
additional rate. The two outermost vectors are defined by the two smallest dot 
products between each of the vectors and the resultant in the V 2 -V ; , plane. 

See the functions between points C and 1) on the flow chart in Appendix A for 
details. The sum of the /i terms do not result in zero torque in the V 2 -V 3 

plane; in fact, the /i terms required to compensate for the torque remove most 

of the , and therefore almost defeat the distribution scheme. For redistri- 
Ai 

Ixi t ion, however, only relative movements are important and, after the summa- 
tion over j is done, the final gimbnl rate additions are 


• • 



V 



( 20 ) 


;j. Throughout the report it is assumed that the final actual and the commanded 
gimbol rates are equal and both are called <’ . ( or /*.) . 
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The disturbance torque caused by the d ^ 


terms are 


T' ■ - ^ ( H oo di d ) 
1 C3*/! V' i l rAi' 


( 30 ) 


and 


t g :i ■ «V A 1> 


(31) 


The inner gimbnl rates also cause a disturbance torque in the V 2 -Vj 


plane, 
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Only the remaining torque needs to lie produced iiy the li terms: 


rp _ p rpt 

CM2 C2 02" 02 


( 3-1) 


and 


T * T _ T’ T" 
CM3 C3 03 03 


(35) 


11 


or 


T 

CM2 





c Vai 


iMj o j) | 


( 30) 


and 


CM3 


T + 
C3 


lilH t ( ou j 


8 H. H 


I' Ai 


♦ HU 




(37) 


We ow apply the pseudo-inverse to net a unique solution. 
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The total outer gimbal rate commands are then 




( 41 ) 
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VI. PROPORTIONAL GIMBAL RATE LIMITING 


Thu gimbal torqucrs will have a definite torque limit, T, .... In order 

UM 

not to exceed this limit, we have to establish the torque demand on the torquers 
due to the gimbal rales. Thi CMG torque is therefore resolved Into the outer 
gimbul coordinate system ( the W, axis is parallel to the outer gimbal axis, the 
W2 axis is parallel to the liner gimbal axis) , 
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where now T is the '.niter gimbul torque (after passing through the inner 

gimbal bear Inn) , T is the inner gimbal torque (which then passes through the 

outer gimbul bearing) , and T ^ Is carried by both the inner and outer gimbal 
bearings. 

Assuming the same torque limit for the iiu.cr and outer gimbal torquers, 
we get for the gimbal rate limits due to the torque limit, 


“tLIM = "tUM* T UM /( "|“'| ) 


(4:i) 


There are also fixed rate limits (<» , ... , /i ,,,) due to other hardware 

1.1M 1.1M 

limits (gimlml rate tachometer limit, vcitage limits, etc.). 
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To reduce the magnitude of the actual torque only, but keep the same 
diruction as the commanded torque, a proportional scaling of all timbal rate* 
by dividing by D is done, where: 


D - MAX 1 
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VII. PERFORMANCE DISCUSSION 


The steering law performance was checked out and verified with a 
hybrid computer simulation. No logie flaws or unanticipated behavior was 
found. 


As implied by the discussion on the desirable gimhnl angle distribution, 
there is no need for a slriet adherence to the ideal distribution. For the outer 
gimbal angles, this means that the time constant ( l/K^) can bo chosen such 

that, while a close adherence for a large V 2 -V : , momentum resultant is 
achieved, small circles of the V 2 -V ;( resultant uboul the V| axis do not 
result in relatively large distributive outer rates. 
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The difference between an outer gimbal angle and its reference eanno. 
exceed ir/2. The maximum distribution rate /3 is then K .»/! , as shown 
in equations ( 28) and ( 29) . 

A more Involved gimbal rate limiting procedure can l>e implemented if 
the dlslriljution rates tend to interfere with the control torque generation. 
Details are shown in Appendix U. 


VIII. CONCLUSIONS 


Parallel mounting of stoploss 1X5 CMG's allows the angular momentum 
steering problem to bo split into a linear and a planar part which simplifies the 
steering law. As a consequence, the steering law can accommodate anv numl>er 
of DG CMG's within reasonable software requirements. The fact that any num- 
ber of CMG's can be accommodated makes treatment of CMC failures a built-in 
property. The capability for clustering any number of CMG's allows the satis- 
fying of any system angular momentum requirement with a few standard sizes. 

The maximum number of DG CMG's necessary in a cluster would be five if the 
next larger size is twice the smaller one (assuming — before failure — no less 
than three CMG's are considered). The stuering-law/purallol-mountod-DG 
CMG combination can therefore be considered a "CMG KIT" applicable to any 
space vehicle where the need for IXi CMG's has been established. 

The logic flow charts of Appendix A are Indiig redesigned and a .'JO per cent 
reduction in logic is anticipated. 
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APPENDIX A 

STEERING LAW LOGIC FLOW CHARTS 
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APPEND I X B 


SUPPLEMENTAL CALCULATIONS 

J, Zero Torque Pisturixmee on V t -Connxme nt. After the summations 
over J were done, the ii nor nimlxtl rale addition (for distribution) was shown 
in equations ( 25) and ( 20') as 
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Yj( ii «» t ) - 2 j( H i ° ,r | ,, i ) 
i i 


(B-l) 


i.e. , no V ( torque is generated for any eombinnlion of II. and <>. values. 


:» l 


I 
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2. Convergence of a to The angle a closely approximates the 

inner reference angle o of equation ( 18) . This can be shown by assuming 

small deviations Ao from o s 

l K 

a ■ or ♦ Aor j ( B-2) 


and 


8 < U R + 11 


( B-:i) 


For A(»j « 1 , wo get 


Sn ''ll — so ' H + ur ''(ll Air ) 
i It ^ i lt^ v i V 


(B-4) 


or 


2(H.A«) = o 
l 


( B-5) 


The angle a becomes 


a - >: 


I llj c(<f R ♦ Anr i)( <r H + L * 11 i L ^ n H + * 


( B-0) 


or 


Of 


kj"i - ”"n * “« [“’r p'i - “"n A, 'ij| 




( B-7) 
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For « l , we get 



(B-W) 

Since V H, An, 0, the difference be tween <> and a., is also zero. 

1 l H 

3. Expanded Gimbal Mate Limiting . 11 the torque demands are large 
when compared with the gimbal rate limits, a more involved limiting procedure 
can lie implemented. When the system is pushed, the outer gimbal rates due to 
the distribution could subtract from the ones needed for the torque generation. 
The inner distribution rates are always small - after the initial transient has 
died out - since the torque generation itself tends to keep the inner gimbal 
angles equal and therefore on the distribution. 

The inner rates are calculated as shown Ijy equation (27) and then 
divided by 


D = MAX ( 1 , l<r j I /<’ ^ . Iif 2 l/« L ♦ ••• • I <** n I ^ (r j ) • 


where 


it 


MIN(,: TUM • " UM> 


( B-‘J) 


This effects limiting, if necessary. 

The outer gimbal rates for the control torque generation arc calculated 
next, using equation ( UH) , 


£ c = IB’I + 


C2 


C:t 


T' 

G2 


- T 


t;:i 


and then divided by 


D - MAX ( 1 , D , D ) 

U (J 


( B- 10) 


where 


D (J - MAX(I/J cl l/» L . I0 C2 I//« L 


l ' i Cn |/ V 


(B-ll) 


and 


0, = MIN (/i 


TLIM 


ii UM ) 


( B- 12) 


This effects limiting, if necessary. Next, the outer distribution rates are 
calculated (equation 29) with 



V 

Li 

j 



To eliminate a torque disturbance, these rates have to be further modified to 


(i A = {E- |B , | t IB'II (£ a 


( B- l:j) 


where K is the n * n identity matrix. 


The i, M 


are then divided by 


I). = MAX ( 1 
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( B- 14) 
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where 



V 1 


D 


H 




( B- lb) 


is the left over capability after the torque generation demand is met. If 

is less or equal to zero, all /i are set to zero. 
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